Femtosecond nonlinear ultrasonics in gold probed with ultrashort surface plasmons 



Vasily V. Temnov 1 *, Christoph Klieber 1 , Keith A. Nelson 1 , Tim Thomay 2 , Vanessa Knittel 2 , 

Alfred Leitenstorfer 2 , Denys Makarov 3 , Manfred Albrecht 3 and Rudolf Bratschitsch 3 
1 Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 
2 Department of Physics and Center for Applied Photonics, 
University of Konstanz, D-784 57 Konstanz, Germany and 
2 Institute of Physics, Chemnitz University of Technology, 09107 Chemnitz, Germany 

(Dated: July 31, 2012) 



o 

(N 

7— I 
> 

X 



£Sp Fundamental interactions induced by lattice vibrations on ultrafast time scales [l]-[3| become increas- 
y— ( ' ingly important for modern nanoscience and technology [4-8]. Experimental access to the physical 
properties of acoustic phonons in the THz frequency range and over the entire Brillouin zone is cru- 
cial for understanding electric and thermal transport in solids and their compounds [9]. Here, we 
i— h \ report on the generation and nonlinear propagation of giant (1 percent) acoustic strain pulses in hy- 
I ^ ■ brid gold/cobalt bilayer structures probed with ultrafast surface plasmon interferometry [8L lid]. This 
new technique allows for unambiguous characterization of arbitrary ultrafast acoustic transients. The 
giant acoustic pulses experience substantial nonlinear reshaping already after a propagation distance 
of 100 nm in a crystalline gold layer. Excellent agreement with the Korteveg-de Vries model points 

i 1 to future quantitative nonlinear femtosecond THz-ultrasonics at the nano-scale in metals at room 

^ ■ temperature. 

■ i-h , Investigations of short wavelength THz acoustic phonons in the nanometer range usually require sophisticated 
^ • methods with high spatial resolution such as X-rays, electrons, or far ultraviolet radiation. However, in metals, 
O . measurements at visible or infrared wavelengths intrinsically bear the potential of nanometer spatial resolution due 
to the characteristic penetration depth (skin-depth) of optical fields of approximately 10 nm. Time-resolved optical 
experiments in metal and semiconductor thin films [l[ provide a detailed picture of the generation of picosecond 
acoustic pulses [llj. However, the detection capabilities in picosecond ultrasonics are generally limited by the largely 
^-j" unknown photo-elastic coefficients leading to complicated dynamics of strain-induced surface reflectivity and ren- 
. dering the characterization of acoustic pulses extremely challenging 0, . Some of these difficulties have been 
i circumvented by using a more straightforward beam deflection technique [13 1 or equivalent interferometric detection 

of transient surface deformations [3J. Nevertheless, the separation of geometrical and photo-elastic contributions to 
H 1 the optical phase still represents a problem. An even more complicated experimental geometry utilizing joint phase- 
^ ' and polarization-sensitive reflectivity measurements at oblique light incidence was theoretically proposed to solve this 
C^ - . problem [14] . Inspired by the performance of the first acousto-plasmonic pump-probe experiments in Kretschmann 
geometry in 1988 [HI and the recent progress in lateral nanofabrication techniques in hybrid multilayer structures 
combined with quantitative nanoplasmonic measurements below the optical skin depth [p, [l6| we have developed fem- 
tosecond surface plasmon interferometry [lOJ . The phase of surface plasmons propagating along the metal-dielectric 
interface is not sensitive to transient surface displacements, but is only affected by the strain-induced photo-elastic 
contribution Q . As demonstrated below, this fact allows us to fully characterize ultrashort acoustic pulses in gold by 
monitoring transient modulations of the surface plasmon wave vector. 

Beyond these important metrological issues, recent studies have revealed the possibility to generate ultrashort 
acoustic pulses by magnetostriction |5[ , to manipulate the magnetization in ferromagnetic semiconductors 0, U3] and 
metals [l8[ by ultrashort strain pulses, and to use acoustic pulses to even generate THz radiation These findings 
render the field of THz acoustics truly interdisciplinary. Pioneering investigations of nonlinear acoustic dynamics in 
dielectric crystals at low temperature 0, Q have revealed clear signatures of shock-front decay into a sequence of 
acoustic solitons after a propagation distance of hundreds of microns. Surprisingly, acoustic nonlinearities have never 
been investigated in metals, although they are used as sources of large-amplitude acoustic pulses for these experiments. 

Here, we demonstrate large-amplitude acoustic strain pulses in a gold generated by femtosecond laser excitation in a 
hybrid gold/cobalt bilayer structure. We perform quantitative measurements by plasmonic pump- probe interferometry 
with femtosecond time resolution [Toj . Figure la shows the experimental configuration of this technique adopted to 
study acousto-plasmonic effects. In contrast to most previous measurements performed with a single polycrystalline 
metal layer on a dielectric substrate, we use a Au/Co/sapphire structure with well-defined crystallographic orientation. 
A hybrid acousto-plasmonic 120 nm gold/35 nm cobalt/sapphire multilayer structure is manufactured by magnetron 
sputtering of a (111) gold layer on top of a hep-cobalt film deposited on a (0001) sapphire substrate. The cobalt 
layer excited by an intense ultrashort laser pulse through the sapphire substrate serves as an efficient and ultrafast 
opto-acoustic transducer. Due to strong electron-phonon relaxation of the time scale of 200-300 fs T9J and the short 
electron mean free path l e ~ 1 nm [2tj| in ferromagnetic materials, hot electron diffusion is much less efficient than in 
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FIG. 1: Acousto-plasmonics in a hybrid gold-cobalt multilayer structure, a, Schematic drawing of the acousto- 
plasmonic pump-probe experiment: surface plasmons propagating at the gold-air interface probe the reflection of acoustic 
pulses generated by the laser-heated cobalt transducer, b, Measured acousto-plasmonic pump-probe interferogram showing a 
pronounced shift of the interference fringes upon reflection of an ultrashort acoustic pulse. The inset illustrates the geometry of 
the plasmonic slit-groove interferometer (see Ref. @,[lp 9 .16] for details), c, Dynamics of the surface dielectric function extracted 
from the plasmonic interferogram. 



noble metals [211) . Therefore, the heat penetration depth only slightly exceeds the skin depth for pump light (typically 
by ~ 50% [Til), see Fig. la. The laser-heated cobalt layer thermally expands and generates an ultrashort acoustic 
strain pulse propagating in both into the gold layer and into the sapphire substrate. Excellent acoustic impedance 
matching between the three layers suppresses acoustic reflection at both the gold/cobalt and cobalt/sapphire interfaces 
(about 10%). As a consequence, the initial shape of a unipolar acoustic pulse follows the spatial profile of deposited 
heat in cobalt. The nonlinear dynamics of ultrashort acoustic pulses in sapphire has been addressed in detail before 
[^LUt- Here, we focus on nonlinear acousto-plasmonic effects in a metal. 

The compressional acoustic pulse in gold r](z, t) = (rii(z, t)— n®)/n® creates a layer of higher ion density rii(z, t) > 
which slowly moves at the sound velocity c s =3.45 km/s in gold in the (111) direction. Since the stationary charge 
separation between electrons and ions in a metal occurs only within the Debye radius rDobye ~ 10 -3 nm, the spatial 
profile of electron density n e (z,t) exactly follows the ionic one: n e (z,t) = rii(z,t). At a photon energy of 1.55 eV 
of the probe pulses the dielectric function e m = e + ie = —24.8 + l.5i in gold is dominated by the free-carrier 
contribution with e ~ — ojp/oj 2 oc — n e . An ultrashort acoustic strain pulse creates a time-dependent spatial profile 

of the dielectric function e (z, t) — e (1 + r)(z, i)) inside the metal, which modulates the surface plasmon wavevector 
k sp = &o\Am/(l + e m): when the strain pulse arrives within the surface plasmon skin depth <5 s ki n = 13 nm at the 
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FIG. 2: Acoustic strain pulses at the gold-air interface, a, The measured strain pulse (blue solid line) launched with a 
peak fluence of 7 mj/cm is in excellent agreement with the theoretical calculation (red line). The inset shows the two terms 
in Eq. (3). b, Simulation of the strain pulses generated in cobalt (blue solid line) and after propagation through a 120 nm gold 
layer based on a linear model including dispersion (dashed line), as well as dispersion and surface roughness (red solid line). 
The black arrow indicates the propagation direction of the acoustic pulses (compressional acoustic strain pulses have positive 
sign throughout the paper). 



gold-air interface: 



5k sp (t) = k ° / r)(z,t)exp(-\z\/6 skln )sgn(z - c s t)dz . (1) 



2|e m |<5. 



skin 



The transformation of the incident compressive pulse into a tensile pulse upon acoustic reflection is described by 
the sign function sgn(z — c s t) in Eq. (1). Figure lb shows the experimental time-resolved plasmonic pump-probe 
interferogram, which depicts the pump-induced changes of optical properties at gold-air interface captured by time- 
delayed surface plasmon probe pulses and recorded in optical transmission through the slit. Pump pulses excite 
the cobalt layer at zero pump-probe delay time. The most striking feature of this pump-probe interferogram is the 
pronounced modulation of plasmonic interference fringes upon reflection of an acoustic pulse at gold-air interface 
at pump-probe delay time of t =39 ps. The dynamics of surface dielectric function in Fig. lc (reconstructed from 
pump-probe interferogram in Fig. lb, see Ref. [1J| for details) shows that the acoustic modulation is due to the strong 
change of the real part of the dielectric function Se (i). The monotonously increasing background signal in both Se (t) 
and Se (t) describes the slow temperature rise at gold-air interface on a time scale of a few tens of picoseconds caused 



by diffusion of hot electrons in gold, an effect extensively discussed earlier [21 1. 

Information of the dielectric permittivity on an ultrafast time scale lets us now reconstruct the acoustic strain 
pulses at the gold-air interface. Assuming that the shape of the acoustic pulse does not change during the acoustic 
reflection, i.e. rj(z,t) = rj(t — z/c s ), we obtain a more convenient expression for Se (t) ~ 2|e m | 2 (5fc sp (i)/fco: 

Se (*) = iiiL / r](t )exp(-|< - t |/T skin )sgn(i - t )dt (2) 

Tskin J — oo 



with 77(f) = rj(z = 0,t — to) and T s km = <5skin/c s . Taking the time derivative de jdt leads to the following Fredholm 
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integral equation of the second kind for the strain rj(t): 

V(t) = W^^T-^— T ^>xp(-|* - t'\/r akin )dt' . (3) 
z\€ m \ at zr skin j_ oa 

Whereas a solution of this integral equation for an arbitrary pulse is a difficult task, it may be simplified for an 
ultrashort acoustic pulse: as long as the duration of the acoustic pulse is shorter than the acoustic propagation time 
through the optical skin depth T s k in = 3.8 ps, the first rapidly varying term in Eq. (3) dominates and the numerical 
solution of Fredholm equation converges after a few iterations. 

In Fig. 2a we show the acoustic strain pulse obtained from the experimental data in Fig. lc by solving Eq. (3). The 
inset in Fig. 2a illustrates the contributions from the two terms in Eq. (3). The first term given by the time-derivative 
de' j dt determines the overall shape of the acoustic pulse. The shape is unaffected by small errors in the determination 
of the slowly varying background described by the second integral term. 

The retrieved acoustic pulse shape may be well approximated by simple linear modeling of acoustic propagation, 
as explained in detail in Fig. 2b. The initial strain pulse generated in the cobalt transducer is modeled with a simple 
exponential exp(— |i|/ro) with tq — 2.5 ps, followed by a weak antisymmetric tail reflected from cobalt-sapphire 
interface (acoustic reflectivity R=-0.11). Upon arrival at the gold-air interface, i.e. after propagation through 120 nm 
gold, the acoustic strain pulse suffers from dispersion [22j . High-frequency components are left behind the pulse 
maximum, as observed in the high frequency ripples in Fig. 2b. Moreover, the acoustic strain pulses experience 
reflection at slightly different times due to nano-scale surface roughness (SR) at the gold-air interface. We used 
atomic force microscopy to quantify SR on all investigated samples. A typical surface roughness of ^2 nm (RMS) 
introduces a distribution of acoustic delay times within 600 fs (FWHM) and smears out the high-frequency ripples of 
the acoustic pulse (Fig. 2b). 

An excellent quantitative agreement between the experimental data and the linear model (Fig. 2a) taking into 
account acoustic dispersion and surface roughness is achieved by setting the heat penetration depth in cobalt (5h oa t = 
Cs Co ^tq to 15 nm with ci Co ^ — 6.3 km/s being the the speed of sound in (0001) direction in hep-cobalt [23j . This 
value is ~1.5 times higher than the optical skin depth for 400 nm pump pulses in cobalt of 9.6 nm, in agreement with 
previously reported results for nickel films [ll[ . 

The measured amplitude of the acoustic pulse ~ 2 x 10~ 3 agrees well with the strain due to thermal expansion of the 
laser-heated cobalt layer, estimated from the amount of absorbed pump energy. Possible contributions due to ultrafast 
magnetostriction [j| would lead to much smaller strain values below 1.5 x 10 -4 in cobalt (24|. It is important to note 
that the acoustic pulses entering the gold layer preserve their duration and exponential shape ~ exp(— t/ro), but due 
to the 1.8 times smaller speed of sound in gold they become higher in strain and spatially contracted. Therefore, 
particularly high strain amplitudes in gold should be feasible, in contrast to sapphire, where acoustic pulses become 
weaker for the same reason(ci s<w '" r6 ' ) = 11.2 km/s in (0001) direction). 

In order to be able to work in the range of relatively high pump fluences we had to use (for technical reasons) pump 
pulses at optical wavelength of 800 nm, which were discriminated from probe pulses at same wavelength by perpen- 
dicular polarization. Using this arrangement we have performed a series of pump-power dependent measurements 
with higher excitation fluence in two samples with an identical thin cobalt transducer of 35 nm, but two different 
thicknesses of the gold layer, 120 nm and 220 nm. We clearly observe nonlinear reshaping of the acoustic pulses in 
both samples (Fig. 3). Larger pulses are moving faster and the reshaping being more pronounced in thicker gold 
film. Both observations indicate the importance of nonlinear acoustic propagation in gold. Our measurements are in 
excellent quantitative agreement with the solutions of the non-linear Korteveg-de Vries (KdV) equation 



di] drj d 3 ij C3 drj 
dt s dz ^ dz a 2pccP dz 



^+-rU + ^r^=o. (4) 



These dynamics are governed by the interplay between the acoustic broadening due to phonon dispersion ui(q) — 
c s q — jq 3 with 7 = 7.41 x 10~ 18 m 3 /s and self-steepening and re shap ing due to the elastic nonlinearity C3 = —2.63 x 
10 12 kg/ms 2 . See Ref. @ for a definition of C3 and Ref. [ID, |26| for linear and higher order elastic constants in 
gold, p — 19.2 g/cm 3 is gold density. Surface roughness is again found to smear out possible high-frequency acoustic 
transients on a sub-picosecond time scale (compare KdV and KdV+SR in Fig. 3). Agreement between the non- linear 
theory and the experiment is simply achieved by setting a single fit parameter, the initial heat penetration depth in 
cobalt <5hcat = 20 nm (corresponding to t$ = 3.2 ps initial pulse duration). This value is again 50 percent larger than 
the optical skin depth in cobalt at 800 nm pump wavelength of 13 nm. 

The quantitative agreement between theory and experiment for 120 nm and 220 nm thin gold layers clearly demon- 
strates that nonlinear reshaping is dominated by acoustic nonlinearities in gold. A possible dependence of the initial 
strain shape in cobalt on pump power obviously plays a minor role. The amplitude of the acoustic strain reaches 



5 



g 120 nm Au, theory 
0.02 




-5 5 

Time (ps) 

q 220 nm Au, theory 




I3 120 nm Au, experiment 
0.02 

0.015 

0.01 

0.005 



-5 5 

Time (ps) 

(j 220 nm Au, experiment 




0.02 



0.015 



0.01 



J \ 20 mJ/cm 2 




■5 5 

Time (ps) 



5 5 

Time (ps) 



FIG. 3: Nonlinear reshaping of giant acoustic strain pulses as a function of excitation fluence in 120 nm and 
220 nm gold films. a,c Numerical solutions of Korteveg-de Vries equation with and without surface roughness (SR). b,d 
Experimental curves showing nonlinear steepening of giant (up to 1 percent) strain pulses, which becomes more pronounced 
for high excitation fluences and for thicker gold layers. The quantitative agreement with the KdV solutions is obtained by 
adjusting a single fit parameter, namely the 20 nm heat penetration depth in cobalt at 800 nm pump wavelength. No fit 
parameters were used to derive the experimental strain pulses. 



peak values r] max = 0.01 percent, corresponding to a compressional stress dp = (? s pr\ m ax = 2.3 GPa before acoustic 
reflection and negative -2.3 GPa tensile stress afterwards. At even higher pump powers, strain pulses up to 1.5 percent 
were obtained but irreversible degradation of the samples was observed, consistent with elevating the temperature of 
the cobalt layer close to the melting point of 1768 K. Control measurements at lower pump fluences confirmed that 
the nonlinear reshaping presented in Fig. 3 is fully reversible. 

Intrinsic damping of longitudinal phonons in gold caused by anharmonic phonon-phonon interactions [22| becomes 
increasingly important for frequencies exceeding 1 THz. However, in this frequency range and for our experimental 
geometry the effect of nano-scale surface roughness entirely masks possible contributions due to phonon attenuation. 
Ultrafast acoustic pump-probe experiments in hybrid multilayer structures with a thicker gold layer and improved 
surface quality [27j may be performed to explore the mechanisms of THz phonon attenuation in metals. 

In conclusion, we have observed nonlinear reshaping of giant ultrashort acoustic pulses in gold at room temperature 
by applying femtosecond time-resolved surface plasmon interferometry. The ability to generate and quantify high- 
strain ultrashort acoustic pulses in hybrid metal/ferromagnet multilayer structures is most relevant for the emerging 
field of ultrafast magneto-acoustics j(| [l7l. [l8l . where the strongest effects were obtained in the high-stain regime 
(Vmax ~ 0.5%) using 200 nm nickel films [l8|~where the non-linear propagation effects are likely to be important. We 
envision using GPa acoustic pulses to switch magnetization in complex magnetostrictive materials such as Terfenol-D. 
This material is characterized by energy barriers between metastable magnetization directions in the MPa range and 
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possesses an exceptionally high magnetostrictive coefficient of ~ 10 3 , which is about two orders of magnitude higher 
than in Ni, Fe or Co 

Theoretical estimates by Abrikosov show that our ultrashort GPa strain pulses can change the distribution function 
of electrons caused by phase-matched electron-phonon coupling and resulting in the trapping of a small sub-ensemble 
of the electrons in the strain-induced deformation potential oc E-p crrn ir] ~ 50 mcV 29]. This effect may stimulate novel 
applications of time-resolved photoelectron spectroscopy leading to the observation of transient strongly correlated 
electron-phonon states in metals on ultrafast time scales. 
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